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The variable solar magnetic activity known as the 11-year solar cycle has the longest history of solar observations. These cycles dramatically affect conditions in the heliosphere and the Earth's space environment. Our current understanding of the physical processes that make up global 
solar dynamics and the dynamo that generates the magnetic fields is sketchy, resulting in unrealistic descriptions in theoretical and numerical models of the solar cycles. The absence of long-term observations of solar interior dynamics and photospheric magnetic fields hinders 
development of accurate dynamo models and their calibration. In such situations, mathematical data assimilation methods provide an optimal approach for combining the available observational data and their uncertainties with theoretical models in order to estimate the state of the solar 
dynamo and predict future cycles. In this presentation, we will discuss the implementation and performance of an Ensemble Kalman Filter data assimilation method based on the Parker migratory dynamo model, complemented by the equation of magnetic helicity conservation and long- 
term sunspot data series. This approach has allowed us to reproduce the general properties of solar cycles and has already demonstrated a good predictive capability for the current cycle, 24. We will discuss further development of this approach, which includes a more sophisticated 


dynamo model, synoptic magnetogram data, and employs the DART Data Assimilation Research Testbed. 


Components of the Solar Cycle Prediction Approach 


Preliminary Analysis of Prediction Solar Cycle 25 Uncertainties 


Variations of solar activity are a result of complicated dynamo processes in the convection zone. We consider this 
phenomenon in the context of sunspot number variations, for which we have detailed observational data during the 
past 23 solar cycles. However, despite the known general properties of solar cycles, a reliable forecast of the 11- 
year sunspot number is still a problem. The main reasons for these forecasting uncertainties are imperfect dynamo 


To characterize the level of solar activity during systematic observations of sunspots, the relative sunspot 
number was introduced by Wolf (1850): R = k(10g + n), where k is a correction factor, 
observing conditions, g is the number of sunspot groups, and n is the number of individual sunspots. 


Previous experience and tests have shown 


depending on based on the dynamo model and sunspot n 


that the EnKF procedure 
umber measurements has 


good predictive capabilities for estimating future solar activity in a time 


range from 7 - 8 years to a whole solar cycle. However, our attempts to 
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a). Panel b) shows the ratio of these data sets. 


Data Assimilation Research Testbed 


Early estimation of properties of Solar Cycle 25 for the 
: sunspot number version 2.0. Panel a) shows two predictions 
: for Cycle 25: 1) prediction obtained for observations that 
include the sunspot number data up to the solar minimum in 
2008 (green curve); 2) prediction obtained using all currently 
available observations up to 2015 (red curve). Blue curve 
shows the best EnkF estimates of the previous cycles based 
on the dynamo model (4) and all available sunspot 


200. 
oof il 
I 
fey 
v WISTS TICs 
years growth time, years 
Sunspot cycle properties: a) asymmetry of the mean shape 
of solar cycles (dotted curve corresponds to Cycle 14, solid 
thick curve to Cycle 19, and thin line to Cycle 23), and b) 
relationship between the cycle growth time and the sunspot 


http://www.image.ucar.edu/DAReS/DART/index.php 
Anderson et al., 2009; Anderson & Collins 2007 
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=> another b observation (red circles). Panel b) shows the model errors of 
SS cycle? number maximum (Waldmeier's rule). The data are from the toroidal magnetic field variations for the case using the 


original sunspot number set. 
Ensemble Data Assimilation 


observations made up to 2008. Panel c) shows these errors 
for the case when all available data (up to 2015) are used. 
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Comparison of the annual sunspot number series: the original (marked ‘original’ solid curve), which 
was used for the solar cycle 24 prediction by (Kitiashvili & Kosovichev 2008), and the ‘revised’ 
version (v2.0, Source: WDC-SILSO, Royal Observatory of Belgium, Brussels) by applying a new 
calibration to the historical and modern data sets (dotted curve, Clette et al. (2015) is shown in panel 
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ne cycle often fail due to 


accumulation of errors. In this section, | consider in more detail the 
possibility of early forecasts of sunspot cycles and present an initial 
estimate for Cycle 25. 
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Testing the prediction capability for solar cycles 16-23. The green curves show the model reference solution. The blue 
curves show the best estimate of the sunspot number using the observational data (empty circles) and the model, for the 
previous cycles. The black curves show the model solution according to the initial conditions of the last measurement. The 
red curves show the prediction results. 


Magnetic helicity shows significantly 
better correlation with future sunspot 
numbers, indicating that the magnetic 
helicity substantially decreases prior to. 
weak sunspot cycles. 
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solar cycle 
Prediction of solar cycles is one of most interesting problems closely linked to dynamo processes inside the Sun. Numerous earlier attempts to predict future solar cycles were mostly based 
on empirical relations derived from observations of previous cycles and provided a wide range of predicted strengths and durations of the cycles. The difficulty is due to our incomplete 
understanding of the physical mechanisms of the solar dynamo and also due to observational limitations that result in significant uncertainties in the initial conditions and model parameters. 
We have developed a relatively simple non-linear mean-field dynamo model, which nevertheless can describe essential general properties of the cycles and the observed sunspot number 
series (such as Waldmeier's rule). Combined with the data assimilation approach, this model provides reasonable estimates for the strength of the following solar cycles. In particular, the 
prediction of Cycle 24 calculated and published in 2008 is holding quite well so far. 


The initial prediction of Cycle 25 shows that this cycle will start in about 2019 - 2020, reach the maximum in 2023 - 2024, and the mean sunspot number at the maximum will be ~90 (for the 
v2.0 sunspot number series) with an estimated error of ~15%. The test simulation runs for early cycle predictions have been performed and identified the following criteria for a successful 
prediction capability: a) the model errors relative to observations should be less than 20% for the last 10 years, and b) the prediction should be performed starting from a period when either 
the toroidal or poloidal field is dominant relative to the other. This corresponds to two moments of time: the polar field reversals shortly after the solar maxima (strong toroidal field and weak 
poloidal field) and during the solar minima (strongest poloidal and weak toroidal fields). 


The next steps of this study are (1) transition from the sunspot number to physical quantities such as synoptic magnetograms and magnetic helicity and (2) using more advanced dynamo 
models that provide a better physical basis for development of a high-fidelity pipeline for solar activity predictions. 
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Simulated test predictions of Cycle 23 using the v2.0 annual sunspot number series. 
Panels a-c and g-i show Cycle 23 estimations (green curves) for the last observing times 
indicated in the figure panels. Black curves show the initial periodic solution obtained from 


the dynamo equations; red circles show the annual sunspot 
show the best EnKF estimate of the model variations. Panels 


number. The blue curves 
d-f and j-I show toroidal 


magnetic field errors of the model for each data assimilation case. 
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